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Abstract 

We present a theoretical study of a spin field-effect transistor realized in a quantum well 
formed in a p-doped ferromagnetic-semiconductor- nonmagnetic-semiconductor-ferromagnetic- 
semiconductor hybrid structure. Based on an envelope-function approach for the hole bands in 
the various regions of the transistor, we derive the complete theory of coherent transport through 
the device, which includes both heavy- and light-hole subbands, proper modeling of the mode 
matching at interfaces, integration over injection angles, Rashba spin precession, interference ef- 
fects due to multiple reflections, and gate-voltage dependences. Numerical results for the device 
current as a function of externally tunable parameters are in excellent agreement with approximate 
analytical formulae. 

PACS numbers: 85.75.Hh, 72.25.-b, 73.23.Ad 



2 



I. INTRODUCTION 



Spintronics has attracted great interest in the scientific community,— advocating the use 
of the spin degree of freedom in electronic devices. Combining this idea with mesoscopic 
transport has stimulated investigations of coherent spin-dependent phenomena. Many pro- 
posed device setups exploit the effect of spin-orbit coupling on the carrier motion . ^I'^i^i^i'^'i^i^ 
The most popular proposal for a coherent spintronic device is the spin field-effect tran- 
sistor (spin FET) proposed by Datta and Das .^°i^^ It consists of a two-dimensional (2D) 
electron gas confined in a semiconductor heterostructure that is attached to two ferromag- 
netic contacts acting as source and drain. Majority-spin electrons injected from the source 
experience a spin precession due to the Rashba effec t ^^i^^i^^ if the magnetization direction 
in the source contact is parallel to the direction of current flow or perpendicular to the 
plane of the 2D electron gas. Tunability of the spin-orbit coupling strength by gate volt- 
ages enables external control of this spin precession and, hence, manipulation of the current 
transmitted at the second ferromagnetic contact. Besides gate-voltage control of the Rashba 
spin-orbit coupling experienced by 2D electrons, which has been successfully demonstrated 
experimentally,— ii^^ efficient injection of spin-polarized electrons from ferromagnetic con- 
tacts into the nonmagnetic part of the spin FET is a key ingredient for device operation. The 
obvious challenges involved in the fabrication of hybrid systems consisting of metallic and 
semiconducting parts, as well as a physical limitation^^ to the amount of spin injection that 
can be achieved in the absence of tunnel barriers at the interfaces, have so far prevented the 
realization of any spin FET device. A possible solution to circumvent these difficulties may 
be provided by the use of diluted magnetic semiconductorsiS' as source and drain. This mo- 
tivates our present study where we investigate transport through 2D hybrid structures with 
ferromagnetic contacts realized in semiconductor heterostructures. An important aspect 
of our work deals with the fact that ferromagnetic (III,Mn)V compounds are intrinsically 
p-doped, implying that currents are carried by holes rather than electrons. The spin prop- 
erties of carrier states in the intrinsically p-like valence bands of III-V semiconductors are 
very different from that in the s-like conduction band. To begin with, several valence bands 
with different effective masses exist. More importantly, however, spin-orbit coupling in the 
valence bands has a more complicated structure than that of conduction-band electrons . 

The aim of this article is a detailed study of an all-semiconductor spin field-effect tran- 
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FIG. 1: Schematic illustration of the proposed device. The two-dimensional hole gas (2DHG) in 
the GaAs part is attached to spin-polarized source and drain contacts, formed by 2DHGs in the 
MnGaAs parts. The gate electrode on the top controls both the carrier concentration and the 
Rashba spin-orbit coupling strength. 

sistor in which the conducting channel is provided by a 2D hole gas (2DHG). The device 
structure we propose is depicted in Fig. ^ A MnGaAs/GaAs/MnGaAs heterostructure is 
overgrow n^^i^^ in the z direction with AIGaAs such that a 2DHG forms at the interface. 
In fact, Mn doping is only required within the quantum well formed at the interface to 
the AIGaAs layer, but the Mn ions outside the well do not disturb. In the proposed setup, 
source and drain are defined by 2D quantum wells accommodating spin polarized holes. The 
carriers in the entire 2DHG are subject to the Rashba effect which leads to spin precession. 
The strength of the Rashba spin-orbit coupling can be tuned by a gate voltage applied to 
the top of the sample.— 

The key ingredient for the functionality of the device is the tunable Rashba spin pre- 
cession in the 2DHG. In Ref. 0, we studied this part within a simplified model allowing 
for an analytic treatment, which enabled us to discuss features that are universal for both 
electron and hole transport. In the present paper, we aim at a more complete numerical 
treatment of transport through the entire device. This includes describing the semiconduc- 
tor valence bands by a 4 x 4 Kohn-Luttiger Hamiltonian instead of restricting ourselves to 
heavy-hole bands only. We take into account the matching properties of the modes at the 
interfaces to the source and drain contacts, which automatically includes interference effects 
due to multiple reflection. As in Ref. we allow for all possible injection angles instead 
of restricting to a quasi-one-dimensional setup. In our analysis we study the purely ballistic 
regime. Scattering due to impurities or to the lateral finite size of the device due to the 
functional form of the spin-orbit coupling Hamiltonian for holes affects the total angular 
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momentum j and not only the spin. This leads to extra damping of spin-precession and 
should be minimized to improve device functionality. Finally, for addressing the response 
of transport to the gate voltage, we employ a capacitive model that takes the variation of 
both the carrier density and the Rashba coupling strength into account. 

The article is organized as follows. We introduce, in Sec. IHl the Kohn-Luttinger Hamil- 
tonian for the ferromagnetic and nonmagnetic parts of the semiconductor quantum well. 
In the following Section IIIIl we describe the mode-matching technique used to calculate 
transmission coefficients for transport through the structure. Our results from numerical 
simulations of transport are presented in Section IIVI After discussing the case of only one 
interface between a ferromagnetic and a nonmagnetic 2DHG, we turn our attention to the 
full transistor geometry with two ferromagnetic contacts. In the latter case, interference 
effects appear due to multiple reflections at the interfaces. In Section El we describe an 
analytical model that approximates the numerical simulations very well and helps us to 
understand how the precession length Lso depends on the Fermi energy E-p in the 2DHG. 
In Section EH we address the response to external gate voltages. The possibility to control 
the spin precession by a single gate voltage Vg that simultaneously modifies Fermi energy 
and Rashba spin-orbit coupling is discussed in detail. 



II. ENVELOPE-FUNCTION DESCRIPTION OF 2D VALENCE-BAND STATES 

In this section, we obtain effective Hamiltonians that describe the valence bands in the 
different regions of the spin FET, namely, the ferromagnetic source and drain contacts doped 
with Mn2+ ions and the undoped nonmagnetic channel in between. In all these regions, holes 
are confined within a 2D quantum well. We use an envelope-function descriptio n^^i^^i^^ of 
the 2D system. The Hamiltonian for the nonmagnetic semiconductor is the sum of a 2D 
quantum-well Hamiltonian if2D plus a Rashba term i^rs, which arises due to the asymmetry 
of the confinement potential Vcon{z)- On the other hand, the total Hamiltonian for the 
ferromagnetic contacts iff is given by plus the term ifpd, which takes into account the 
coupling between the p-like valence holes and the half-filled d-shell Mn^+ ions with spin 
S = 5/2. 

Bulk systems host heavy- and light-hole bands with total angular momentum j = 3/2. 
These bands are degenerate at the band edge and are well separated, due to spin-orbit 
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coupling, from the split-off bands with total angular momentum j = 1/2. In quantum wells, 
each of these bands is transformed in a sequence of quasi-2D subbands, and the degeneracy 
between heavy and light-hole bands is lifted. In the following, we are interested in the 
situation at low carrier concentrations such that only the lowest subbands, one heavy-hole 
band (HHl) and one light-hole band (LHl), have to be taken into account. This is the 
simplest realistic case that occurs when the triangular quantum well is narrow enough to 
sufficiently lift the energy of higher subbands. Inclusion of higher energy subbands, in 
particular HH2, would be straightforward, but would only lead to higher order corrections. 
Indeed, the low hole density typically present in experiments implies that the only occupied 
propagating modes are in HHl, whose shape is influenced by the band mixing2& with LHl 
and in a negligible way with HH2. Coupling with HH2 would affect the shape of LHl, and 
therefore only evanescent modes. In the basis of total angular momentum, 

|l) = |j = 3/2,j, = 3/2) 

|2) = |j = 3/2,j, = -l/2) 

|3) = |j = 3/2,j, = l/2) 

|4) = |j = 3/2,j, = -3/2) 



the Hamiltonian if2D reads^^ 



2D 



f hh d \ 
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\ d* hh J 



(2) 



with 



hh = E^^'^ + 

2mhh|| 

Ih = E't^i + (3) 

d = -^^Nk'^e-'^" . 
Here we have adopted the momentum-space representation in polar coordinates for the wave 
vector k|| = {k cos a, k sin a) in the 2D plane. The quantities E^^^ and E^^^ in Eq. (jS)) are 
the subband-bottom energies deriving from the solution of the Schrodinger problem for the 
triangular well, and the factor takes into account the scalar product between the envelope 
functions Ji^i^i'i^i for the HHl and LHl subbands. We observe that in Eq. ^ there is no 
coupling between LH and HH subbands with the same sign of jz since the corresponding 
matrix element is proportional to the vanishing integral {f^^^\kz\f^^^)- 
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FIG. 2: Dispersion of the lowest quasi-2D subbands (HHl and LHl) of a semiconductor quantum 
well when no exchange field and no spin-orbit coupling are present. The subband-bottom energies 
are computed using eEz equal to 4 x 10^ eV/m and vertical masses of rrihhz = 0.38mo and mih^ = 
0.09mo. 
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The effective masses appearing in Eq. are given by 

^hhll = m/ (7i + 7) 
mih|j = m/ (71 - 7) 

where the two coefficients 71 and 7 = (72 + 73)/2 are the Luttinger parameters,— taken 
in the so-called axial approximation.— In Fig. |21 we show the HHl and the LHl subband 
dispersion relations as a function of the magnitude k of 2D wave vector. 

The Rashba spin-orbit coupling arises from the structural inversion asymmetrjsi^, due to 
an asymmetric confining potential. It is described by the Hamiltonian i/rs = /5(k x E) ■ j. 
Here, E is the electric field due to the confining potential K:on(-2), and /? is a material 
parameter ."2^^ In our system, E = i?^ z, and H^-s reads 



/ 








\/3p— ia 
2 ^ 





\ 












V3 - 
2 ^ 


-ia 




\/3 pia 
2 












\ 





\/3 pia 
2 








I 



(5) 



In Fig. 121 we show the HH and the LH subbands for the nonmagnetic semiconductor as a 
function of fc, obtained by diagonalization of = H21) + -f^rs- We see that the splitting 
at small wave vectors is linear for the light-hole subbands but cubic for the heavy-hole 
subbands. 



The Hamiltonian for the ferromagnetic-semiconductor part is given by Hf = H2D + 
Hj-s + Hpd- We use a phenomeno logical description of the ferromagnetic semiconductor, 
in which local moments with S = 5/2 from Mn^"*" ions are antiferromagnetically coupled 
to the itinerant holes .^^i^^i^^ In a mean- field treatment, combined with a virtual- crystal 
approximation , ^"^i^^i^^i^^ itinerant holes experience an exchange field h = </pd^Mn(S), where 
the average Mn-ion spin S has direction fi = (cos sin 6, sin sin 6, cos 6), Nun is the doping 
concentration, and Jpd describes the coupling strength. This exchange-coupling field is 
accounted for in i/pd = h ■ a, where the spin matrices a in the basis given in Eq. are 
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FIG. 3: The first subbands (HHl and LHl) of a semiconductor quantum well in the presence of 
Rashba spin-orbit coupling, which removes the spin degeneracy. The coupling constant {(3Ez) is 
0.15 eV nm. 
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FIG. 4: The ferromagnetic-semiconductor 2D quantum-well subbands for the case of magnetization 
perpendicular to the 2DHG plane. The interaction constant is Jpd = 0.06 eV nm^, the Manganese 
concentration is A^Mn = 1 nm~^, and {(3Ez) = 0.15 eV nm. 
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In Fig. 13] we show the HH and the LH subbands in the ferromagnetic-semiconductor contacts 
as a function of fc, obtained by diagonalization of H^. The splitting of the two heavy-hole 
subbands for small values of k leads to full polarization at low densities. The magnetization 
direction in the Figure is perpendicular to the 2DHG plane, n = (0, 0, 1). 



III. QUANTUM STATES FOR HOLES PROPAGATING THROUGH THE SPIN 
FET 



We calculate coherent transport through the spin FET using the scattering formalism 
described, e.g., in Ref. |38l It relates the current to transmission amplitudes for scattering 
states defined in the contacts. To obtain these, proper matching of wave functions at inter- 
faces is required, that we describe in this section. The imposed conditions at the interface 
are the continuity of the wave function and conservation of the component of probability 
current that is perpendicular to the interface. To illustrate the subtleties associated with the 
second condition, let us consider an interface between a ferromagnetic and a nonmagnetic 
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region at x = xq. The continuity equation for the current is 



(9) 



with 




which derives from the operator relation = ^[H,x]. Note that, due to the presence of 



spin-orbit couphng, the derivative of wave functions needs not to be continuous at x = Xq. 
Instead, Eq. Q guarantees current conservation. 

In both the ferromagnetic and nonmagnetic regions, four different channels are available, 
z = 1, . . . , 4, associated with the four-dimensional Hilbert space of the valence-band subspace 
under consideration. Let us consider a wave incoming from the ferromagnetic region with 
wave vector = k] (cos a, sin a), which is in the ith subband. The wave is partially reflected 
at the interface to the nonmagnetic region. The wave function in the ferromagnetic source 
electrode is then 



where Xn the eigenf unctions of H{ in k-space representation, the normalization factor 
f^(k„) is the velocity expectation value computed for the state vector Xni^n), and rj„ are 
the reflection coefficients to be determined from the matching. The wave vectors k^ = 
(cosaj^, sina^) of the reflected waves are determined by the following three conditions: 

a) the modulus k^ is the solution of the implicit equation e^(fcj^) = Ep, where el^{k) is the 
dispersion relation; 

b) the angles are derived from the continuity of the momentum parallel to the interface 
due to translational invariance along that spatial direction; 

c) among these solutions we allow those that satisfy {3fJ[u^(k^)] > & 53[f^(k^)] = 0} or 
{$5[f^(k^)] > 0}. The two possibilities correspond to propagating and evanescent modes, 
respectively. 

Similarly, the wave function in the nonmagnetic region is 




(11) 




(12) 
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where ti^n are the transmission coefficients, and all the other quantities in Eq. p2p have the 
same definitions given for the corresponding quantities in Eq. The eight coefficients 

Ti^n and ti^n are determined by the two conditions of continuity of the wave function and 
current conservation, Eq. Q, given that ip^'^ has four components in the total momentum 
space. 

We emphasize that it is important to include all modes even when some of them are 
evanescent since transmission and reflection at the interface are influenced by tunneling into 
classically forbidden channels. 

IV. RESULTS OF NUMERICAL SIMULATIONS 

For our numerical simulations we take Jpd = 0.06 eV nm^, A^Mn = 1 nm~'^ and S = 5/2. 
The Rashba term is characterized by {(3Ez) = 0.05 eV nm, and the Fermi energy is taken 
as Ep = 0.08 eV, which corresponds to a hole density of 4 x 10^^ m~^. We assume that the 
leads' magnetization direction is either perpendicular {6 = 0) or within the plane (6' = 7r/2) 
of the quantum well. 

To clarify the underlying physics, we will approach the full spin FET design step by step. 
First, we consider Rashba spin precession for holes transmitted through a single interface 
between a ferromagnetic and a nonmagnetic 2DHG for a flxed injection angle. Then, we 
include a second interface with the magnetization of the source and drain electrode being 
parallel, keeping the angle of incidence for spin-polarized holes still flxed. As a result of the 
Rashba effect, the total transmission will oscillate as a function of the channel length in the 
nonmagnetic 2DHG. Finally, we take the full 2D nature of the device into account by adding 
up the current contributions for all injection angles. 

Let us start by considering a single interface between a ferromagnetic and a nonmagnetic 
2DHG. In Fig. El we show the spin-up and spin-down currents as a function of the distance 
from the interface in the case of magnetization perpendicular to the plane (n = (0, 0, 1)) and 
perpendicular injection {a = 0) of spin-up current for two different values of Ep. We flnd 
that both the spin-up and spin-down current density oscillates with modulation length Lso, 
indicating Rashba spin precession in the nonmagnetic region. With increasing Fermi energy, 
the oscillation length decreases. This result is in clear contrast to the case of the spin FET 
based on electrons, where the spin precession length is independent of the Fermi energy.^^ 
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FIG. 5: Spin components of the current density in a nonmagnetic 2DHG plotted as a function of 
the distance L from the interface with a ferromagnetic 2DHG. Results are shown for a half-metallic 
ferromagnetic contact having magnetization direction n = (0,0,1). The Fermi energy is equal to 
0.09 eV in a) and to 0.045 eV in panel b). It is apparent that the period of current oscillations is 
proportional to 'i/E-p. For both spin directions, the current is normalized to the incident hole flux. 

Moreover, the modulation length -C/so(«) depends also on the injection angle a (not shown 
in Fig. [SJ. For a realistic sample, integration over all possible injection angles is required. 
We will see later that, after integration, the overall modulation length is given by that for 
perpendicular injection, -Lso(O). 

In Fig. ini we show results for the case of magnetization direction in the ferromagnetic 
2DHG being n = (1/^2, l/v/2, 0). For in-plane magnetization in the half-metallic contact, 
the amplitude of the current modulation in the nonmagnetic 2DHG depends on its azimuthal 
angle (f). The largest oscillation amplitude occurs for = 0, i.e., when the magnetization di- 
rection is perpendicular to the interface. No oscillations exist for = 7r/2, because majority 
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FIG. 6: Same as Fig. El a) but for the case where the magnetization direction in the ferromagnetic 
2DHG (and, hence, the quantization axis for spin components of the current) is equal to fi = 
(l/-v/2, 0), i.e., hes in the 2DHG plane. Note the diminished amplitude of current oscillations 

which would disappear altogether for n = (0, 1,0). 

spins injected into the nonmagnetic 2DHG are then eigenstates of ifrg. 

We now turn to the simulation of the spin FET transistor, consisting of a finite strip of 
nonmagnetic 2DHG with two interfaces with ferromagnetic contacts, one at a; = and the 
other one at x = L. Now we have to apply the mode-matching procedure for each interface. 
In Fig. [71 we plot the total transmission through the entire device as a function of the width 
L of the nonmagnetic region. We find a modulation of the transmission with modulation 
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FIG. 7: The transmission probability of spin-up electrons for the case of two interfaces separated by 
a distance L. The high-frequency oscillations are due to resonances arising from multiple reflections 
between the two interfaces. It turns out (see below) that such features tend to be smeared out 
when the transmission is averaged over the injection angle. Parameters are the same as in Fig. 
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FIG. 8: The total current density as a function of the quantity L/L^q in the 2D system treatment, 
where L is the channel length and L^o is the total modulation length. The magnetization in the 
contacts is perpendicular to the plane of the 2DHG. 

length Lso, which is due to Rashba spin precession. This modulation is superimposed by 
fast oscillations of the order of twice the Fermi wavelength Ap which are due to interference 
effects from multiple reflection within the double-barrier structure. '^^ As we will see below, 
these fast oscillations will be almost always smeared out after integration over the injection 
angles a, i.e., they will not appear in real 2D devices. Only in the limit of very low hole 
densities, remnants of these oscillations will be visible. 

Finally, we take into account the full 2D nature of the device, i.e., we add up the current 
contributions for all injection angles. We assume an isotropic angular distribution of injected 
holes since all our simulations are performed in the linear response regime. As a result, the 
transmitted current density J / Jq is given by the formula 



where T„,„(«) is the transmission probability from channel n to m for holes injected at 
an angle a. Since the modulation length Lso(a) of the transmission is a-dependent, one 
might expect that the integration washes out the effects of the spin precession. It turns 
out, however, that oscillations are still visible, although damped. The modulation length of 
the resulting oscillations coincides with that for perpendicular incidence, -Lso(O). In Fig. |S1 
we show the result of the integration operation corresponding to Eq. (fT^ for magnetization 
direction in the contact 2DHGs being fi = (0, 0, 1). 




(13) 
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V. ANALYTICAL RESULTS FOR TRANSPORT 



In this section we show that most features of the numerical results presented in the previ- 
ous section can be understood within a simplified analytical model. In particular, we derive 
analytic expressions for the precession length Lso(a) and the total current density. A similar 

was used to discuss universal fea- 
tures of hole and electron spin precession. The approximate formulae derived in this section 
clearly show how the precession length depends on system parameters and, therefore, allow 
for a deeper understanding of the underlying physics than looking at the purely numerical 
results presented in the previous Section can provide. 

The model is developed following some approximations that are justified a priori by 
physical considerations and a posteriori by the comparison between analytical and numerical 
results. First, we make use of the fact that for typical parameters only the lowest heavy- 
hole subband is occupied. As an approximation we can, therefore, omit all non-conducting 
subbands from our model. Furthermore, we assume perfect transmission at the interfaces, 
i.e., we neglect reflection. All Hamiltonians are now represented as 2 x 2 matrices, using the 
basis (H)) restricted to vectors |1) and |4). The off-diagonal matrix elements are obtained 
from perturbation theory for the degenerate case.— The nonmagnetic 2DHG region is then 
described by 

/i^p / 1 \ , . . [ e-*3° 
H, = 7r:—\ +^{f3^E^)k'\ , (14) 



2mhhi| I 1 / I e*^" 

where /3h is proportional to the spin-orbit coupling of holes and is different from the (3 defined 
in Eq. ®. 

The corresponding eigenenergies are 



e,,2{k) = -^e±{P^E,)e (15) 
2mhh|| 



with eigenvectors 



^^''^ V2 (^^ie^3" ' ' ^^^^ 
The spin splitting of the eigenvalues ()15|). together with the conservation of the wave vector 
parallel to the interface, implies the presence of a double refraction phenomenon^Mi where 
a hole wave incident on the interface from the ferromagnet gives rise to two transmitted 
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waves in the nonmagnetic 2DHG having shghtly different wave vectors. Their magnitudes 
ki^2 are obtained from the imphcit equation ei,2{k) = Ep. Typically, spin-orbit coupling can 
be treated as a perturbation, which means that we can linearize the expression of ki 2 in 
the spin-orbit coupling strength, and arrive at /ci 2 = A;o =F AA;/2. Here kQ is the Fermi wave 
vector in the absence of spin-orbit coupling, and 

/ 2mhh|| \ ^ 



Ak 



(17) 



\ J 

which explicitly depends on the Fermi energy. The corresponding angles of the transmitted 
waves' propagation direction with the interface normal are found, again in the limit of weak 
SO coupling, to be 

«i,2 = tto i (A^/2A;o) tanao ; (18) 



where ao is defined by 

fcp sin a = ko sin ao . 
Hence, the transmitted hole is described by the wave function 



^^^_^^iki{xcosai+ysmai) _|_ ^^^^^1/02(2: cos «2+y sin 02) 



(19) 



(20) 



By assuming a perfectly transparent interface, we can compute the coefficients Ci 2 simply by 
matching the wave functions in the ferromagnet and in the nonmagnetic semiconductor. At 
the interface at x = L to the second ferromagnet, for the case of its magnetization pointing 
in positive z direction, only the |+) component will be transmitted. Hence, the outgoing 
state in the right ferromagnet reads e^'^pC^cosa+s/sma) cos[AA;L/(2 cosa;o)]|+). As a result, the 
transmission probability is 

To,<p{a) = cos^ — — , (21) 
[cosaoj 

where we have used the relation AkL/2 = 7, and the dependence on a is through ao 
via Eq. (fT^ . In a similar way we can obtain the transmission probabilities for arbitrary 
magnetization direction in the ferromagnetic 2DHGs. The transmission probability for in- 
plane magnetization reads 



T 



TT/2,(f> 



cos 



7 



cosao 



+ sin^ [Sao — 0] sin^ 



7 



cosao 



Finally, we can write the transmission for arbitrary magnetization direction as 

Te,ct>{a) = cos^ 6^0,^ + sin^ 6^ T^/2,</, . 



(22) 



(23) 
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Equations. ()21II23|) cease to be valid once one of the transmitted states in the nonmagnetic 
2DHG becomes evanescent, i.e., is totally reflected. This condition defines critical angles 
ttc,{i,2}, that in the limit of weak SO coupling read ac,{i,2} = ^o/^f =F |A/c/A;f ~ k^/k-p = ac- 
We note that very similar formulae for the transmission can be obtained for electrons.-^ 
From Eqs. (j?T|l and (f^^ we find the precession length 

^ , , 27!- / \^ cos an /^.N 

{(JhEz) V2"^hh||/ 

which depends on both the injection angle a and the Fermi energy Ep. The physical reason 
for the latter dependence is the cubic spin splitting of the heavy-hole subband. For a 2D 
device under consideration in this paper, one has to integrate over all injection angles, see 
Eq. (fT^. The final result reads then 

J/ Jo = I cos^ ^F(7) + sin^ 6 

X sin^ + F(7) cos^ + ^(7) cos(20) } , (25) 

where Jq is the injected current density, and the functions -^(7) and G(7) are defined as 

Fh) = U% ^os a cos^^) da, (26) 
^(7) = i /i cos a sin^ (3a) sin^ {^) da . (27) 



A good analytical approximation for -^(7) and G{'y) is given in Ref.|2^ The result Eq. ^ 
describes damped oscillations of the current density as a function of the length of the non- 
magnetic part, where the modulation length is given, in agreement with our numerical 
findings, by the precession length for perpendicular injection, Lso(a = 0). 



VI. GATE- VOLTAGE MANIPULATION OF THE CURRENT: CAPACITANCE 
MODEL 

In this section we consider the response of the hole spin FET to external gate voltages 
and point out important differences in its behavior compared to the electron version. With 
standard densities of about 10^^ m~^, only the lowest spin-split 2D HH subbands are oc- 
cupied. Their k^ spin splitting leads to an inversely linear dependence on the Fermi energy 
Ep for the hole precession length Lg^. Hence, variation of gate voltages will modify L^^ by 
changing, at the same time, the asymmetry of the hole confinement in the 2DHG and the 
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Ca Cbg 
2DHG 



Cq 



X 



BG 



J_ Vd=V3 



FIG. 9: Linearized capacitance model for the influence of top and back-gate voltages on the spin 
FET. The 2DHG is located at point A, Vq is the voltage applied at the top gate, Vbg that of the 
back gate. The voltages at the drain (Vd) and source (Vs) contacts are assumed to be equal since 
we consider the linear response regime. 

Fermi energy. We analyze how the device performance changes when, instead of only a top 
gate voltage, both top gate and back gate voltages are applied. 

We model the effect of the top and back-gate voltages through a linearized capacitance 
model as shown in Fig. El The capacitance per unit area between the top gate and the 
point A, where the 2DHG is located, is Ca = Co^GaAs/c^A ~ 2 x 10~^ F/m^ for an effective 
distance = 50 nm between top gate and 2DHG, while the capacitance associated to the 
variation of the back-gate voltage is Cbg = Ca/2 for a distance of 100 nm from the 2DHG. 
The capacitance between the 2DHG and the source and the drain is well approximated by 
the quantum capacitance Cq = mhh||e^/7r^^ ~ 8 x 10^^ F/m^ and is due to the finite density 
of states in the 2DHG.— The variation of the voltage at the point A and of the electric field 
Ez read then 

Ca , Cbg 



dVA 

dE.. 



Ca + Cq + Cbg 
Cq + Cbo_^^^^ 



dVc,+ 



Ca + Cq + Cbg 
dV^G) 

eo^GaAs 

(Cq + CBG)\CAdVG - (Ca + CQ)dVBG] 



dV, 



BG 



A 



eoeGaAs(CA + Cq + Cbg) 

while the variation of the Fermi energy reads 

-e {CAdVG + CBGdVBc) 



dE-p = —edVA 



Cq + Cbg + Ca 



(28) 



It is clear that variation of only Vq (i.e., keeping dVBc = 0) simultaneously changes E^ 
and Ep. In order to leave the Fermi level pinned, we have to manipulate both top and 
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back-gate voltages such that c/Vbg = — (C'a/C'bg) c^Vg- Results of our simulation for the 
case that only the top-gate voltage is varied {(IVbg = 0) is shown in Fig. ^1 Here the 
current density obtained from Eq. (|13p is plotted as function of the gate voltage Vg for the 
case when the magnetization direction in the contacts is perpendicular to the 2DHG. Spin- 
precession-induced current oscillations are clearly visible. The oscillation period actually 
varies with changing Vg due to the induced variation of Ep, as mentioned in Sec. IIIII Note 
that, for the parameters of Fig. [TU[ the effect of the gate voltage on the position of subband 
bottoms is negligible since the electric field is changed only by a few percent from its initial 
value Ez ~ ep/eoCGaAs- (Here p is the hole density.) For comparison, we show results for 
the case where both top and back-gate voltages are varied simultaneously such that the hole 
density in the 2DHG remains unchanged in Fig. ^2 Here the precession length changes only 
due to the gate-voltage-induced variation of the structural inversion asymmetry, measured 
here by the electric field E^. Current oscillations have then a larger period as function of 
Vg than in the case where Vbg is kept constant. In both cases, however, a clear modulation 
of the current as a function of gate voltage Vg is obtained. These are slightly damped due 
to the superposition of current amplitudes for all possible angles of incidence. 

0.85 
° 0.75 
0.65 
0.55 

-0.1 -0.2 -0.3 -0.4 -0.5 
Vg(V) 

FIG. 10: The modulation of the current density through the spin FET when the top-gate voltage 
is varied. The magnetization direction in the ferromagnetic 2DHGs is given by n = (0,0, 1). The 
Fermi energy for Vg = is Ep = 0.09 eV, the Rashba coupUng term {(3Ez) = 0.05 eV nm. 
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0.55 

-0.2 -0.4 -0.6 -0.8 

Vg(V) 

FIG. 11: The modulation of the current density in the spin FET when the top-gate voltage and 
the bulk-gate voltage are simultaneously varied in a way such that E-p is unchanged. The contact 
magnetization points along fi = (0, 0, 1). 

VII. CONCLUSIONS 

We have performed careful numerical and analytical studies of transport through a p-type 
all-semiconductor spin FET. The design of such a device would overcome problems associ- 
ated with the fabrication of hybrid devices involving metal-semiconductor contacts. Despite 
the more complicated nature of spin splitting in 2D valence-band states, clear current mod- 
ulation as a function of device parameters such as the width of the nonmagnetic region are 
observed. Using a phenomenological model for the action of external gate voltages, we have 
shown the possibility of current manipulation as was envisioned in the original electron spin 
FET proposal by Datta and Das.— Our numerical simulations that were performed for re- 
alistic sample parameters, as well as analytical formulae reported in this work, should serve 
as a useful basis for experimental realization and functional optimization of a p-type spin 
FET. 
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